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Abstract: We present a non-perturbative formalism for measuring defect free en- 
ergies (monopole mass or vortex tension) in three-dimensional SU(2)+adjoint Higgs 
models. Starting from twisted, translation invariant boundary conditions, we per- 
form a change of variables that allows us to express the defect free energies in terms 
of 't Hooft loops. We propose that the defect free energies can be used to distinguish 
between phases in this model, and also more generally in other gauge field theories 
where no local order parameters exist. In the case of monopoles, our construction 
can also be used in four- dimensional pure-gauge SU(2) theory, where it gives the 
monopole mass in the maximally Abelian gauge without the need of actually fixing 
the gauge in the simulation. Moreover, the expression is manifestly independent of 
the choice of the Abelian projection as long as it is compatible with the classical 
't Hooft-Polyakov solution. 
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1. Introduction 

In many cases, phase transitions can be described in terms of topological defects. In 
some simple models, such as the Ising model [0] and super- Yang- Mills theories 0, an 
exact relation is known between this dual formulation and the original formulation of 
the theory. However, even without such an exact mapping, this dual description can 
be very useful, in particular since it gives a way of finding natural order parameters 
for the phase transition. 

If the phase transition is associated with a spontaneous breakdown of a global 
symmetry, essentially any local operator that transforms non-trivially under the sym- 
metry group can be used as an order parameter. However, if the broken symmetry 
is a gauge invariance, this approach does not work, because only gauge-invariant 
operators have non-zero expectation values. Although sometimes, for instance in the 
electroweak theory 0, |^, the transition predicted by perturbation theory is in fact 
only a smooth crossover, often a well-defined transition exists. In these cases, a dual 
picture in terms of topological defects can often be the most natural description. 



There are two possible ways of constructing an order parameter based on topolog- 
ical defects. In theoretical studies, it is convenient to consider the creation operator 
/i of a defect 0, ||, |^, ||, ^ 0, [11|. In the low-temperature phase, its expectation 
value (fi) vanishes, but if the defects condense in the high-temperature phase, it be- 
comes non-zero. Thus it behaves as an ordinary order parameter in a spontaneous 
symmetry breakdown, but with an inverted temperature, and often it can indeed be 
associated with a symmetry of the dual theory that is broken in the high-temperature 
phase. Therefore it is often called a disorder parameter ^j. However, this approach is 
only useful if the defects are line-like and can therefore be interpreted as world-lines 
of particles of the dual theory. 



More generally, one can use the defect free energy |T^, |13| |1^, |1^, |T6[ as an order 
parameter. In the special case of line-like defects, the tension, i.e. free energy per 
unit length, is simply the inverted correlation length of fi, and thus vanishes at the 
transition point. Using the defect free energy as the order parameter also avoids 
the problem that in numerical simulations (/i) always vanishes because of the finite 
lattice size. 

In the Abelian Higgs model, Nielsen- Olesen vortices have been used successfully 
to study the phase transition. In Refs. [^, ||, the disorder parameter and the cor- 
responding U(l) symmetry were discussed, and in Ref. |jl3l the tension of a vortex 
was measured and shown to be an order parameter. The purpose of this paper is 
to extend these studies to SU(2) theories with one or two adjoint Higgs fields by 
constructing explicitly the observables that measure free energies of vortices and 
monopoles. 

Even though there is no spontaneous gauge symmetry breakdown in pure-gauge 
SU(2) theory, it is still believed that the non-perturbative dynamics of the gauge 
fields can give rise to similar effects. Both vortices fl^ and monopoles [Q have been 
suggested as possible explanations for confinement. While the vortex tension can be 
expressed simply as the expectation value of a 't Hooft loop fl^ and the procedure 
of measuring it is well understood |114|, Il5|, [16|, giving a non-perturbative definition 



for the monopole mass is more involved. It requires defining a composite adjoint 



field with an Abelian projection |]T^, after which monopoles can be defined in the 
same way as in the Georgi-Glashow model pp] , pT| . Disorder parameters based on 
monopoles have been studied in Refs. |^, 0, ^, but they depend on the choice of 
the Abelian projection, refer to solutions of classical field equations and don't have 
the same elegance as expressing the vortex tension in terms of a 't Hooft loop. 

In this paper, we present a gauge invariant, non-perturbative formalism for 
studying vortices and monopoles in SU(2) theories with adjoint Higgs fields. In 
Sect. D, we discuss topological defects and defect free energies in general terms. 
Sect. ^ contains the derivation of a gauge-invariant expression for the residual U(l) 
gauge field in the broken phase. In Sect. ||, we discuss the two-Higgs model and 
review the definition of the vortex tension in terms of a 't Hooft loop. This acts as 



an introduction to the definition of tlie monopole mass in tlie Georgi-Glasfiow model, 
wliicli is carried out in Sects. ^. In Sect. |], we sliow liow tlie same construction can 
be used in pure-gauge theories with Abelian projections. 

2. Topological defects 

In continuum, topologically stable non-trivial solutions of the field equations, i.e. topo- 
logical defects, exist when the vacuum manifold of the system has non-trivial topol- 
ogy. They are characterized by a winding number, which is non-zero if a defect is 
present. The mass of such a defect can be defined simply as the energy of the field 
configuration. When the symmetry is restored, the distinction between a defect and 
the vacuum disappears, and we can say that the mass of the defect vanishes. 

However, in quantum systems and in classical systems at finite temperature, the 
state of the system is not given by a single field configuration, but rather by a density 
operator or an ensemble of configurations. In both cases, a useful description for the 
state of the system is given by the partition function, which has the form 

D(f)exp{-S). (2.1) 

Here S = J d^xC is the action and in the quantum case, we have performed the 
Wick rotation. In the thermal case, S = jSH, but for simplicity we assume that the 
temperature is absorbed in the parameters of the Lagrangian C. 

The partition function is a path integral over a large number of field configu- 
rations, almost none of which are solutions of the equations of motions. For each 
configuration, the winding number of the whole system is still a well-defined quan- 
tity. When it is zero, there can be localized objects that at a suitable length scale 
have the characteristics of a topological defect, but the total number of defects and 
anti-defects (those with a negative winding number) must be equal, and therefore 
they can be thought of as defect-antidefect pairs created by thermal fluctuations. 
However, if the total winding is non-zero, the number of defects and anti-defects is 
not equal, which means that the configuration contains, in addition to the thermally 
generated defect-antidefect pairs, true topological defects. Note that in general it 
is not possible to determine which of the defects are the true ones and which are 
members of the thermally generated pairs. 

From the point of view of local observables, it does not make any difference in 
the thermodynamical limit if we use instead of the canonical partition function ( |2.1j ) 
a microcanonical one Zn, where the integration is restricted to configurations with a 
given total winding number n. However, the value of Z„ is different for each n. In 
fact, it is more useful to consider the free energy 

F„ = -lnZ„, (2.2) 



because in the classical limit where the integral is dominated by its saddle point, the 
free energy difference 

AF = Fi - Fo (2.3) 

is exactly the classical mass of the defect. We will call AF the free energy of a defect. 
Note that it is defined without recourse to the saddle-point approximation. 

It is not obvious how to use Eq. (|2.3| ) in lattice field theories, because the topology 
of the field configuration space is very different in a discrete space. For instance, the 
Abelian Higgs model has two different lattice formulations, the compact and the 
non-compact one, and only the non-compact version contains topological defects. 
Thus we will first have to show that in any particular case we are discussing, the 
winding number can really be defined. Even that is not enough, because in lattice 
simulations, it is only possible to measure expectation values 

(C) = Z-^ f D^Oexp{-S), (2.4) 

and not free energies directly. Therefore, we will also have to rewrite Eq. ( |2.3|) in 
terms of expectation values. 

3. Residual U(l) invariance 

Let us start by considering in general an SU(2) gauge theory with an adjoint scalar 
field $, which may be either fundamental or composite. We will now derive an expres- 
sion for the Abelian gauge field corresponding to the residual U(l) gauge invariance 
that remains if the field $ breaks the SU(2) symmetry. A similar construction has 
previously been carried out in continuum [^ and on a lattice in Refs. p3| , p4 . 



In this section, we assume that the system is defined on a three-dimensional 
lattice, which may also be a single time slice of a four- dimensional lattice. The 
gauge field is represented by SU(2) matrices Ui{x) defined on the links {x,x + i) 
between the lattice sites x G {0, . . . , A^ — 1}^, and the adjoint scalar field $ is defined 
on lattice sites. 

The theory is invariant under gauge transformations 

$(f) -^ A"^(f)$(f)A(x), 

Ui{x) -^ A\x)Ui{x)A{x+t), (3.1) 

where A(x) is an SU(2)-valued function defined on the lattice sites. In particular, it 
is always possible to gauge transform $ to the 2;-direction. 

To simplify the discussion, we neglect the configurations in which $ either van- 
ishes or is proportional to o"3 at any lattice site. Because they have a zero measure 
in the partition function, this does not change the results, and our final results work 
well even in these special cases. This allows us to define the unit vector 

$ = $($2)-i/2_ (3 2) 



We start the construction by making a gauge transformation that diagonahzes 
$, i.e. turns $ into a^ at every point. That is accomphshed by 

R{x)(xi{a^ + ^{x)). (3.3) 

Here and in the following we use the notation ex to show an equality that is true 
up to a real and positive factor. Our final expressions will be independent of these 
factors. In Eq. ( |3.3|) , the normahzation is chosen in such a way that R is an SU(2) 
matrix. 

We can now define the transformed gauge field 

tJi{x) = R\x)Ui{x)R{x + i). (3.4) 

Since the corresponding value of $ is proportional to a^, -R is nothing but a gauge 
transformation into the unitary gauge. However, this does not fix the gauge com- 
pletely, because of the residual U(l) invariance. More precisely, a gauge transforma- 
tion with a matrix A induces a transformation 

Ui{x) -^ A'^{x)Ui{x)A{x + i), (3.5) 

where 

A = i?Ui?A, (3.6) 

and 

RA(xt{a^ + A^^A). (3.7) 

The new transformation A is unitary, 

A^A = R^A^RR^ARa = 1, (3.8) 

and its determinant is one. It is also diagonal, since it leaves a^ invariant by con- 
struction, and therefore it must be of the form A = exp(zAcr^). This is the residual 
Abelian gauge transformation. 

Let us then show that the corresponding gauge field is given by the phase of the 
diagonal elements of t/j. We define the projection operators 

and the projected fields 

V,,±{x) = Pjj,{x)P±. (3.10) 

The field 1^+ {Vi^J) corresponds to the upper left (lower right) component of the 
total gauge field f/j. Furthermore, if we define 

n± = i(i±<l), (3.11) 



we can write Vi^±{x) = K^{x)Vi^±{x)R{x + i), where 

Vi,±{x) = n±(f)f/,(f)n±(f +?). (3.12) 

Under a gauge transformation, the fields V^ -t transform as 

Vi^±{x) -^ A^{x)Vi^±{x)A{x+i) = exp {±i (A(f+?) - A(f))} Vi,±{x), (3.13) 

which means that 



ai = arg Tr Vi^^ 



arg Tr V^ _ 



(3.14) 



indeed behaves hke the Abehan gauge field. Thus we can also define the magnetic 
fiux density as 



a. 



,, ^ argTry,,+(f)\/,-+(f+?)V;|+(f+j)V}|+(f) (3.15) 

= arg Tr n+ (f ) f/, (f ) n+ {x + i) Uj {x + i) n+{x+i+j) UJ (x+j) n+ (x+j) U] (f ) . 

This final version is gauge invariant and contains no reference to the fields in the 
unitary gauge. 

Using Eq. ( |3.15| ), we can define the magnetic charge inside a lattice cube as 



Cm{x) 



-^ Y^ tijk {ajk{x+i) - ajk{x)) . 

ijk 



(3.16) 



This number is an integer and can be non-zero. 



4. 3D two-Higgs model 



Let us first consider a model with two adjoint Higgs fields, $ and x- The Lagrangian 
is 



C 



4 



2H 



ag 



5^ M - ^ReTi Ui{x)Uj{x+i)UJ{x+j)UJ{x] 
+ ^ 2a [Tr <^^{x) - Tr $(x)f/,(f)<l'(f +?)[// (f) 

i 

+ Trx^(a^) - Tr x{x)Ui{x)x{x+i)UJ (x) 
+m|a^Tr (!?\x) + a^A$(Tr ^^ 
+mJa^Tr x^(x) + a^A;^(Tr x^)^ + a^r/Tr $\2 + a^r]2 (Tr <l>xf ■ (4.1) 



If $ is non-zero, the two-Higgs model looks very much like the Abelian Higgs 
model. In Fig. ^ we show the plot of the susceptibility of Tr x"^ at the point where x 
becomes non-zero in perturbation theory (see the Appendix for numerical details). 
The maximum susceptibility does not increase with the volume, and therefore the 
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Figure 1: The susceptibility of Tr^^ near the perturbative transition point. 

field Tr x^ does not behave non-analytically. The same apphes to any local observ- 
able, and therefore we are forced to consider non-local order parameters, and in 
particular the defect free energy (^.31) . Because in this model, the defects are line- 
like vortices, it is convenient to define the tension as the free-energy per unit length 

of a vortex 

AF 

lim . (4.2) 

L—*00 Li 



T 



The construction of this order parameter is well known and dates back to late sev- 
enties [|T^ ^, but we will rederive it here, because the derivation of the monopole 
mass in the one-Higgs model follows similar lines. 



4.1 Winding number 

We start by deriving the lattice expression for the winding number. We will follow 
the corresponding derivation in the case of the Abelian Higgs model ||2^ (see also 
Ref. 10). Since we already know the analogue ( |3.14| ) of the Abelian gauge field, we 
only have to find the corresponding Higgs field. If we define using Eq. 



X 



R^xR, 



(4.3) 



we find that it transforms as 



X -^ A^xA. 



(4.4) 



If we now project to the lower left matrix element of x 

h = P^xP+. (4.5) 

we find the transformation law 

h = P.RxRP+ = RU^x'n.+R -^ i?AA^n_xn+Ai?A = A^hA = exp{2i\)h. (4.6) 

Thus h behaves as a charged scalar field. Note, however, that the charge of h is two, 
and therefore the transformation A = tt leaves it unchanged. This transformation is 
the remaining Z2 symmetry that is left after h gets an expectation value. 

Now, we need a gauge invariant expression for the difference of the Higgs phase 
angle at neighbouring lattice sites 

= TYx{x)U_{x)Ui{x)U^{x+i)x{x + i)U+{x+t)UJ{x)U+{x). (4.7) 

Defining the Higgs phase angle 7 by /i oc exp{ij{x)) and using Eq. ( |3.14| ), we find 
(cf. Eq. (10) in Ref. 



6^ = arg A, = [j{x+i) - -f{x) - 2a,{x)]^ , (4.8) 

where we use the notation [X]„ = X + ^nnx, with such an nx that [X]^ G (— vr, vr]. 
To define the winding number we have to subtract the magnetic flux. This leads 
to the expression 

Yij (x) = 5, (x) + 6 J (x+i) - 6, (f + j) - 6 J (x) - 2aij (x) , (4.9) 

which is obviously a multiple of 27r and thus the winding number riij = {27r)~^Yi,j is 
always an integer. The value of Yij is only defined modulo An, and therefore one can 
only distinguish between even and odd values of riij. If it is odd, there is a vortex 
going through the plaquette. 

The definition of the winding number can be extended to more complicated 
curves than just single plaquettes in a straightforward way by adding all the (5's and 
aj's along the curve, and the result is additive 

nc+c = nc + nc mod 2. (4-10) 

This equation implies that vortices cannot end. 

4.2 Vortex tension 

Now that we have the expression for the winding number, we have to restrict the 
path integral to configurations with a given winding. The most natural way of doing 



it is by boundary conditions. In lattice simulations, the boundary conditions are 
almost always chosen to be periodic, 

$ (f + Nj) = $ (f ) , x{x + Nj) = x{x), Uk (x + Nj) = Uk{x). (4.11) 

This implies periodicity also for 6i and Oj, and therefore the total winding is even, 
i.e. there are no vortices. 

Besides their simplicity, one reason for using the periodic boundary conditions 
is that they minimize the finite-size effects, since they guarantee that the actual 
boundaries will be invisible to the physics. This is particularly important when 
a defect free energy is measured, because if there are physical boundaries on the 
lattice, their contribution to the free energy can dominate over it. However, any 
boundary conditions that preserve the translation invariance of the system avoid 
these problems, and therefore it is enough to have boundary conditions that are 
periodic up to any symmetries of the theory. We can, for instance, modify the 
boundary conditions ( |4.11| ) by 



TT f , AT \ I -U2{x,y,z), if y = l, 



Changing the sign of U2{x, 1, z) does not change the Lagrangian (|4.1j ) and therefore 
there are no boundary effects. Moreover, 6i are still periodic, since they always 
contain f/j as well as t//. However, ay{x, 1, z) changes by vr, and if we calculate the 
winding of the whole xy plane according to Eq. ( [4.9|) , ay{x, l,z) appears once, and 
thus the total winding changes from even to odd. This shows that these non-periodic 
boundary conditions lead to configurations with one vortex in z-direction. The free 
energies in Eq. ( pl3| ) are thus given by 

-\nJvU,V^Vxexp [-X^AhJ , (4.13) 



F„ = -InZ 



and the subscript n E {0, 1} indicates which of the boundary conditions ( |4.11| ), ( [4.12| ) 
are used. 

Next, we have to write T = AF/L in terms of expectation values. Let us 
first notice that when we calculate the action S = ^£2H; the boundary condition 
( [4.12|) only affects the Wilson term. Instead of saying that we are using non-periodic 
boundary conditions, we could say that the boundary conditions are periodic, but 
the action is changed by 

AS = ^J2^TU^iN,l,z)U2{l,l,z)UliN,2,z)uUN,l,z). (4.14) 

This corresponds to a change in the integration variable in the path integral, which 
should be taken into account when measuring any observable, but since we are now 



only interested in the tension, we do not have to worry about that. The translation 
invariance of the boundary conditions implies that we could equally well choose any 
coordinates Xq, Hq and write 

A5=^VTrf/i2(xo,yo,^). (4.15) 

This allows us to write 

r=-hm'"<^''P'-^^». (4.16) 

L— >oo Li 

This form is not, however, particularly suitable for numerical simulations, because 
the observable exp(— AS*) and the integration weight exp(— 5) have only very little 
overlap. One solution would be to measure derivatives of T. For instance, 

dT 



X, 



a' 



iV2((Trx')i-(Trx')o), (4.17) 



where the subscript indicates whether the expectation value is calculated in the one- 
vortex or the zero-vortex ensemble, this can then be integrated to yield T. The 
drawback of this method is that it cannot be used to measure the absolute value of 
T, only its differences. 

Another possibility is to introduce a real number e G [0, 1] and define non- 
physical ensembles with the action 5* + eAS. We can define 

F,^-\n f VUiV<!>Vxexp{-S - eAS), (4.18) 

and it is then easy to write the tension as 

T = -y de^ = lj M^S),, (4.19) 

where (. . .)^ means an expectation value calculated using the action S + eAS*. The 
tension can now be calculated in numerical simulations by measuring the integrand 
of Eq. ( |4.19| ) at a large number of different values of e. Note that although ensembles 
with non-integer e are not physical, the final result T is. This method can be improved 
by using multi-histogram techniques |T^, ^, ^|. An integration scheme in which 



e = 1 but plaquettes are added to Eq. ( [4.15|) one by one has also been suggested 
recently [[T6| . 



The shift ([4.15| ) of the action is a 't Hooft loop |JT^ |25| that is made stable by 



closing it via the periodic boundary conditions. Its relation to the twisted boundary 
conditions ( [4.12|) was first discussed in Refs. [^ 0. 



In pure-gauge SU(2), the vortex tension has been measured in Refs. [0, |T5|, |l| 
Instead of having a closed loop, the sum over z in Eq. ( 4.15|) can also be restricted 
to a shorter interval {zi, . . . ,Zf} 

{fi{zf)fxiz,)) = Z-' J VU{D^Vx exp ( -5 - ^ X] ^^ f/i2(a;o, l/o, z) J . (4.20) 



10 



This technique was used in Ref. |^| to study the monopole-antimonopole interaction 



in the pure-gauge SU(2) theory. If the space is interpreted as 2+1 dimensional and 
the z direction as time, this describes a creation of a vortex at time Zi and its 
annihilation at time Zf [^]. This leads to a dual picture of the theory in which 



the fundamental degrees of freedom are vortices and the phase transition has an 
interpretation of a spontaneous symmetry breakdown of a magnetic symmetry for 
which (/i) defined by 

(^)2 = hm (/i(z)/i(0)) (4.21) 

2— >00 

is an order parameter. 

5. 3D Georgi-Glashow model 

Let us now consider a system with only one adjoint Higgs field. We will use the 
Lagrangian 

+ ^2a [Tr<l>2(f) - Tr <l>(f)t/i(f)$(f + i)f//(f) 

i 

+m2a^Tr$2(x) + a3A(Tr$2)2_ (5.1j 



The topological defects in this system are the 't Hooft-Polyakov monopoles pO, 
|2l|| , which are point-like objects. Thus we can define the mass of a monopole using 
Eq. ( |2.3|) simply as 

M = AF. (5.2) 

Mean field theory predicts that monopoles are massive in the broken phase and 
massless in the symmetric phase. Therefore the mass would be a useful order pa- 
rameter for the phase transition. However, the masslessness of monopoles in the 
three-dimensional theory would not imply that there are infinite correlation lengths 
in the symmetric phase, because a monopole is a point-like object and therefore its 
mass is not related to the correlation length of any operator of the dual theory. In four 
dimensions, creation and annihilation operators of monopoles have been discussed in 



Ref. 1 11 1 using a different approach. 



5.1 Boundary conditions 

Again, our strategy is to find boundary conditions that preserve the translation 
invariance of the system but force the total magnetic charge of the lattice to be one. 
Periodic boundary conditions 

^x + Nj) = <l>(f), Ukix + Nj) = Uk{x), (5.3) 



11 



are obviously ruled out, because they don't allow non-zero total charge. Instead, 
we need more complicated boundary conditions, and we can use as a guidance the 
classical continuum monopole solution 



^(x) 



XkCTk 



r 



Mx) 



^ijkXjOk 

2r2 ■ 



(5.4) 



If we move from one boundary to another, we reverse the sign of the coordinate Xj, 
and the fields transform as 









$ 



-a^-^aj. 



Ai > ajAiCTj. 



This suggests the boundary conditions 



^{x + Nj) 



-aj^{x)aj, 



Uk{x + Nj) = ajUk{x)aj 



(5.5) 



(5.6) 



It is straightforward to see that this is a symmetry of the Lagrangian ( p. 11 ), and 
therefore does not break translation invar iance. 

The boundary conditions ( ^.6|) imply that 
n+(a; + Nj) = ajll-{x)aj and, consequently. 



aij{x + Nk) = —aij{x), 



(5.7) 




i.e. the direction of the magnetic flux is re- 
versed at the boundaries. Thus, when we 
cross the boundary of the lattice, we enter 
a charge-conjugated copy of the same lattice \ 
from the opposite boundary. This is however -^ 
not enough to guarantee that the total charge 
of each of these copies is non-zero. 

To determine the total charge, we cal- 
culate the magnetic flux through half of the 
boundary. The curve defined in Fig. || separates the boundary into two halves, and 
we can calculate 



Figure 2: The curve used in Eq. 

to separate the boundary of the lattice 

into two halves. 



flux 



arg 



Tr n ^ M 



N-l 



[X 



arg Y[ JjTrV,,+ (m + A^i)V;|+(m + M), (5.8) 

ijk cyclic n=0 



where the final form follows from the fact that the V^^+ commute. To calculate this, 
we need to know the boundary conditions induced by Eq. (|5.6|) for R: 



R{x + Nx) = — (Ti-R(x)cri, 

R{x + Ny) = -(T2i?(s)cT2, 

R{x + Nz) oc a3R{x)[^{x),a3]a3. 



(5.9) 
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Then it is easy to calculate the corresponding boundary conditions for Vi^±: 
Vi^+{x + Nx) = aiVi^^{x)ai, 

V,,+ {x + Nz) oc a3[a3, <l>(:r)]K_(x)[<l(f+z), (Jajas. (5.10) 

The case i = 3 does not contribute, because 

1iV^^j^{nz + Nx)Vlj^{nz+Ny) = Ti aiV'i^^{nz)aia2Vl _{nz)a2 

= i:i%_{nz)Vl_{nz) = I, (5.11) 

where we have used the fact that o"iO"2P± = icr^P^ = ±iP±. For the horizontal 
segments we find 

Tr Vi,+{x+Ny)Vl^{x+Nz) 

ocTrVi,_(f)ai[(T3,<l>(f+x)]V;|_(f)[<l>(x),a3]CTi, 

Tr V2,+{^+^^)^l+i^+^^) 

ex Tra2[crs,Hx)]V2,4^)mx+y),a^]a2Vl_ix), (5.12) 

where we have written ni = x for notational simplicity. Because 

(5.13) 

and because Vi^_ projects to the lower right component, we obtain 

Tr Vi^^{nx + Ny)V{j^{nx + Nz) oc exp {—i [9{{n + l)a;) — 9{nx)\} , 

Tr V2,+{ny + Nz)Vl^{ny + A^x) oc exp {i [d{{n + l)y) - d{ny)\} , (5.14) 



where $i + i$2 oc exp(z6'). In Eq. ( |578| ) all other phases cancel except those at the 
corners, and we obtain 



-'5 



flux = argexp {i [6{Ny) - e{Nx)]} = arg(-l) = n, (5.15) 



where we have used the boundary conditions for $ in Eq. ( ^.6| ). 

Thus we have shown that with these boundary conditions, there will be a flux vr 
through each of the halves of the boundary, and since the boundary conditions force 
the direction of the flux to be opposite, this means a total flux of 2tt from the lattice. 
Note that because the flux is only defined modulo 27r, it is only possible to force the 
total magnetic charge to be either even or odd. 
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5.2 Monopole mass 

Using the boundary conditions (|5.6| ), we can now define the effective mass of an 
isolated monopole. For any boundary conditions, the free energy is defined by 



F = -In 



J DU^D<!>expl-J2^Mx)] 



(5.16) 



If we denote by Fq and Fi the free energies for ensembles with boundary conditions 
( [5.3|) and (p. 61), respectively, we can use Eq. (|5.2| ) directly to define the monopole 
mass M. 

In order to write M as an expectation value, we will now transform the one- 
monopole system in such a way that the difference is moved from the boundary 
conditions into a shift in the action. Let us consider a field configuration that satisfies 



the conditions (5^), and apply the following (large) gauge transformation 



A 



1, iix,z<N 
iai, ii X < N, z > N 
ids, if X > A^, z < N 
io"2, ii x,z > N 



(5.17) 



This will change the boundary conditions into 

^x + Nj) = -a2Hx)a2 = <^*{x), Uk{x + Nj) = a2Uk{x)a2 = U;{x) 

everywhere else except at the edges of the lattice, where 

Us{x, N,N-1) = -f/3*(x, 0, AT - 1), 
Ui{N-l,N,z) = -Ui*{N- 1,0, z), 
Ui{N-l,y,N) = -Ui*{N-l,y,0). 



(5.18) 



(5.19) 



The C-periodic boundary conditions ( |5.18|) have been discussed before in Ref. [^ in 
a different context. 

When calculating the free energy ( |5.16|) , we can change the boundary conditions 
everywhere to Eq. ( 5.18| ) by redefining 



Usix, N,N-l)^ -f/3(x, N,N-1), 
Ui{N-l,N,z) ^-Ui{N-l,N,z), 
U,{N -l,y,N) ^ -U,{N -l,y,N). 



(5.20) 



This changes the Wilson term in in the Lagrangian ( ^.1|) . The total effect is to flip 
the sign of the Wilson term at three of the edges, and consequently the one-monopole 
free energy is given by 



Fi = -In / DUiD^ex^p - 



X / 



(5.21) 
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with C-periodic boundary conditions ( 5.1S ). The change of the action is 

AS = —rRelYT:TU23ix,N-l,N-l) 

Af-1 N-l \ 

+ 5^ Tr f/i3(iV - 1, y, iV - 1) + ^ Tr Uu{N -l,N-l,z)]. (5.22) 

y=0 2=0 / 

Because these boundary conditions preserve translation invariance, we could as 
well choose any xq = (xq, yo, Zq) and write 

4 ^^-^ 



AS = -^Re ^ Tr ^723(0;, yo, Zq) 



^ \ x=0 



N-l N-l 

+ J^ Tr f/i3(xo, y, ^o) + X^ Tr Uuixo, yo, z) j . (5.23) 

y=0 2=0 



Note how similar this expression is to Eq. ([4.15|). The only difference is that we 



have three intersecting 't Hooft loops instead of just one. Again, we emphasize that, 
because Eq. ( ^.21|) is equivalent to Eq. (|5.16|) with the translation invariant boundary 



conditions Eq. (|5.6|) , the choice of xq does not affect any observable, and in particular, 
it does not fix the location of the monopole on the lattice. 

Without AS*, Eq. ( |5.21|) would be precisely the ordinary free energy ( |5.16|) with 



C-periodic boundary conditions ( |5.18|) . As they obviously imply 

Vi^+{x + N3) = a2Vi,_(f)a2, (5.24) 

the flux through the curve in Fig. ^vanishes (cf. Eq. (|5.11|) ). Consequently, we can 
instead of the periodic boundary conditions (|5.6| ), use the boundary conditions in 
Eq. (|5.18|) to define Fq. This leads to the expression 

M = -ln(exp(-A5)), (5.25) 

where AS" is given by Eq. (|5.23|) and the boundary conditions by Eq. (|5.18|) . 
To measure M, we can use the derivative with respect to m?, 

^ = T/((Tr$^),-(Tr$^), (5.26) 

where the subscript tells whether the expectation value is calculated in the zero- 
monopole or one-monopole system, and V = a^N^ is the volume of the system. If 
we assume that at large enough m^, i.e. in the symmetric phase, the monopole mass 
vanishes, we only have to measure DM /dm? at sufficiently many values of m? and 
integrate. 
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Like in the case of a vortex, we can also define for e G [0, 1] 



F, = -ln I DUiD^ exp j - ^ Cm{x) - eAS j . (5.27) 

Using Eq. (|5.27|) , we can now write the nionopole mass as 

M = de—^ = / de{AS),, (5.28) 

Jo c'e Jo 

where the subscript e indicates that the expectation value must be measured in the 
presence of the insertion —eAS* as in Eq. ( |5.27| ). This gives us the absolute value of 



M, but with the cost that we have to measure expectation values at non-physical 
values of e. 

6. 4D pure gauge SU(2) 

Because of the large amount of literature on condensation of monopoles as an expla- 
nation of confinement in pure-gauge SU(iV) theory |Q |10|, |1^, ^ |23[, it is interesting 



to see whether our formulation can be applied to this case as well. 
The Lagrangian of the 4D pure gauge theory is 



Ad = /3$^ (l - ^ReTTU,,{x)] . (6.1) 



In order to define monopoles, we will also have to specify the Abelian projection |]T9| 
we are using, i.e. how $ is related to the gauge fields. The original suggestion by 
't Hooft [O] was to choose an operator O that is a product of the link variables and 



to define O = Oq + i^, where Oq is proportional to the unit matrix. 

In order to use our definition for the monopole mass, the boundary conditions 
must be chosen in such a way that the induced boundary condition for $ is compatible 
with Eq. ( |5.6| ). However, for most choices of O this is not possible at all. This is 
obviously the case with all definitions of O that are products of spatial links f/j, 
because then the boundary conditions (p.6|) imply 



0{t,x + Ni) =aiO{t,x)ai, i.e. ^{t,x + Ni) = ai^{t,x)ai, (6.2) 



in contradiction with Eq. (p.6|). As a consequence, these Abehan projections are also 
incompatible with the classical monopole solution ( p. 41 ). 

There are suitable ways of defining $, though, and one of them is the most 



popular choice of an Abelian projection, namely the maximally Abelian gauge |]19, 



[23[| . It means finding the gauge transformation A(t, x) that maximizes the functional 

RuAG = Y, Tr a-^t/^(t, x)a^Ul{t, x), (6.3) 



fl,t,X 
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where Ui{t,x) = A^{t,x)Ui{t,x)A{t,x + i) and Uo{t,x) = A^{t,x)Uo{t,x)A{t + l,x). 
Defining $(t, x) = A(t,x)cr^A"'"(t,x), we can write Eq. ( |6.3|) as 



-Rmag — 7 ^ 



t,x 



Tr |)(t, x)Uo{t, f)l'(t+l, x)t/o^(t, f) 
+ ^ Tr <l(t, f)f/,(t, x)l'(t, x + z)f//(t, 



X) 



(6.4) 



which must be maximized with respect to the field $(t, x). It is obvious that if 
$ maximizes -Rmag in the region x E {0, . . . ,N — 1}^ and U^ and $ satisfy the 
generalization of the boundary conditions ( ^^ , 

l>(t, f + Nj) = -cjj^{t, x)aj, U^{t, X + Nj) = ajU^it, x)aj, (6.5) 

then $ maximizes -Rmag everywhere. Thus the maximally Abelian gauge is compat- 
ible with the boundary conditions (p.6[). 

We will now assume that $ is given either by the maximally Abelian gauge 
condition or some other way that is compatible with Eq. (|6.5| ). The definition of the 
magnetic field can be generalized from Eq. (|3.15|) trivially by simply using the same 
definition at every time slice. Consequently, we can perform the gauge transformation 
( |5.17| ) and the field redefinition (|5.20| ) or equivalently simply take Eq. (|5.25|) and write 

ln(exp(-A^)) ^gg^ 

Lt is 



M 



lim 



where the boundary conditions in the spatial directions are C-periodic ( |5.1^ 
the length of the lattice in the time direction and 



x=0 



(6.7) 



N-l N-1 

+ ^ Tr f/i3(t, xo, y, zq) + J^ Tr Uuit, xq, yo, z) 

y=0 z=0 

As before, the mass can be measured in practice using the multi-histogram tech- 
niques 0, |2^, ^ or plaquette by plaquette [|16| . Note that Eq. ( |6.6| ) does not refer 
to $ at all and therefore, the tedious task of numerically maximizing -Rmag is avoided 



altogether. Moreover, as long as the definition of $ is compatible with Eq. (|6.5| ), the 
mass of a monopole is independent of the precise choice of the Abelian projection. 

Let us now compare our prescription with that used in Ref. 0. There, instead 
of the monopole mass, the disorder parameter (fi) was measured from the correlator 
{ji{ti)jj,{tf)). The correlator was defined by shifting the vertical plaquettes UiQ in the 
action at t = t,- 



Uio iU ,x) = Ui {U , x) f/o (t, x+ i) U} (ti + 1 , x) Ul {U , x) 

-^ U,{ti,x)Uo{ti,x+t)UJ{ti + l,f)f/Ji .(ti + l,x)U^oiti,x) 



{6.t 
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where Uci,i is the gauge field of a static classical Abelian monopole, and correspond- 
ingly at t = tf. By redefining the spatial links Ui, 

f/d,,(t,f)f/i(t,f)-^[/i(t,f), (6.9) 

the shift can be moved to the spatial plaquettes, thus making direct comparison 
possible. Leaving aside the problems caused by the fact that the redefinition ( |6.9| ) 
changes the boundary conditions and with most choices of the Abelian projection 
also the value of $, we can write the resulting shift in the action as 

ASit) = -/5 J] J]Re f^TTU:^{t,x) - ^TTU,,{t,x)\ , (6.10) 



Kj 



where f/; = f/i.[/,. 



The structure of Eqs. (|6.7] ) and ( |6.10|) is similar, and realizing that in the full 



non-perturbative treatment there is no reason for f/ci,j to be a solution of the classical 
field equations, we can interpret Eq. (|6.7|) as simply a particularly convenient choice 
of ?7ci,j- For instance, it makes it much easier to choose suitable boundary conditions 
that preserve the translation invariance and allow a non-zero total magnetic charge. 
The measurement of the monopole mass has also been discussed in Ref. [^ 



There the non-zero total magnetic charge was achieved by using fixed boundary 
conditions. However, this breaks the translation invariance, and boundary effects 
can therefore be significant. 

7. Conclusions 

In this paper, we have presented a way of measuring free energies of topological de- 
fects in SU(2) gauge theories on a lattice. Our approach was to impose "twisted" 
boundary conditions, which preserve the translation invariance and minimize bound- 
ary effects but guarantee that the number of topological defects inside the system is 
non-zero. By changing the integration variables in the partition function, the free 
energy of the defect can then be expressed as an expectation value, which makes it 
possible to measure it in numerical simulations. In the case of the two-Higgs model, 
we derived an expression for the vortex tension in terms of a 't Hooft loop using this 
technique. 

The main result of this paper is the derivation of an analogous expression for 
the mass of a 't Hooft-Polyakov monopole in the Georgi-Glashow model, defined 
as the free energy difference between configurations of total magnetic charge one 
and zero. Again, the construction was based on manifestly translation invariant 
boundary conditions, which guarantees that the boundary effects are minimal. By 
a gauge transformation and a redefinition of the fields, we rewrote the free energy 
as a path integral with C-periodic boundary conditions and with an insertion AS, 
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which consists of three intersecting 't Hooft loops. In addition to avoiding problems 
with boundary conditions, our prescription leads to a simpler final expression for the 
monopole mass than the ones appeared previously in the literature. 

We believe that the monopole mass (|5.25| ) can be used to distinguish between 



the phases of the Georgi-Glashow model, because it is zero in the SU(2) phase and 
non-zero in the U(l) phase. More generally, we conjecture that if the symmetry 
breaking structure of a phase transition in a gauge theory allows topological defects 
and it is possible to construct an observable that measures the free energy of such 
a defect in the same way as Eqs. ( |4.16| ) and ( |5.25| ) measure the free energies of 



vortices and monopoles, the point at which the defect free energy becomes non-zero 
defines a transition point between the phases. Although this does not necessarily 
imply a singularity in the partition function, it suggests that in these cases the phase 
transition is not likely to be a smooth crossover. 

We also showed that the same construction can be used in pure-gauge theory to 
observe condensation of Abelian monopoles. In this framework, it is not necessary to 
fix the gauge or carry out the Abelian projection in the simulation. This simplifies 
the procedure significantly and means that the observable, and therefore also the 
resulting monopole mass, is manifestly independent of the choice of the Abelian pro- 
jection, provided that it is compatible with the classical 't Hooft-Polyakov monopole 
solution. Thus, the ambiguities related to the differences of various Abelian pro- 
jections can be avoided by adopting Eq. ( |6.6| ) as the definition of condensation of 
Abelian monopoles. 
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A. Computational Details of the Two Higgs model calculation 

The configurations for this calculation were generated using a combination of the 



Kennedy-Pendleton [^ heatbath algorithm and the over-relaxed 0, ^ algorithm 
along with a Metropolis accept-reject step for the pure gauge sector and the scalar- 
gauge interactions and a Metropolis step for the scalar interactions. Measurements 
are carried out every 500 compound sweeps. In Table |1], we list the thermalization 
and number of measurements carried out for the lattices used in Fig. 0. 
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V/a" 


Number of thermalization steps 


Number of measurements 


123 


50,000 


200 


163 


50,000 


200 


243 


10,000 


40 



Table 1: Run parameters for lattice sizes presented in Fig. 1. 



As outlined in [^, if we define 

P(i)(f) = aTrx?i), 



(A.l) 



where i indicates the z-th configuration, (the fields computed on the lattice are ^/ax) 
and 



then the susceptibility in Fig. ^ is defined as 

where the notation {(D) represents an average over configurations. 



(A.3) 
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